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ABSTRACT: In this article, we describe a method to increase the density of field-effect transistors in 


the framework of an absorptive single pole single throw switch. In the framework of the approach we 
consider to manufacturing of inverters in heterostructure with specific configuration. Several required 
areas of the heterostructure should be doped by diffusion or ion implantation. After that, the dopant 
and radiation defects should by annealed in the framework of optimized scheme. We obtain 
dependences of optimal value of annealing time on several parameters. We also considered a method 
to reduce the amount of stress mismatch in the considered heterostructure. We consider an analytical 
approach to analyze mass and heat transport in heterostructures during manufacturing of integrated 
circuits with account mismatch-induced stress. 
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INTRODUCTION 


In the present time several actual problems of 
the solid state electronics such as increasing of 
performance, reliability and density of elements 
of integrated circuits: diodes, field-effect and 
bipolar transistors are intensively solving [1-6]. 
To increase the performance of these devices it 
is attracted an interest to determine the 
materials with higher values of charge carrier 
mobility [7-10]. One way to decrease the 
dimensions of the elements of integrated 
circuits is manufacturing them in thin film 
heterostructures [3-5,11]. In this case, it is 
possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12,13] and 
development of epitaxial technology to 
improve these materials (including analysis of 
mismatch induced stress) [14-16]. An 
alternative approach to increase the dimensions 
of integrated circuits are using of laser and 
microwave types of annealing [17-19]. 


In the framework of the paper we introduce an 
approach to manufacture field-effect transistors. 
The approach gives a possibility to decrease 
their dimensions with increasing their density 
framework a absorptive single pole single throw 
switch. We also consider the possibility to 
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decrease mismatch-induced stress to decrease 
the quantity of defects, generated due to the 
stress. In this paper, we _ consider a 
heterostructure, which consist of a substrate 
and an epitaxial layer (see Fig. 1). We also 
consider a buffer layer between the substrate 
and the epitaxial layer. The epitaxial layer 
includes into itself several sections, which were 
manufactured by using other materials. These 
sections have been doped by diffusion or ion 
implantation to manufacture the required types 
of conductivity (p or n). These areas became 
sources, drains and gates (see Fig. 1). After 
doping it is required to annealing of dopant 
and/or radiation defects. 


Fig. 1a: Structure of the considered switch [15] 
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Epitaxial Layer 


Fig 1b: Heterostructure with a substrate, epitaxial layers and buffer layer 
(view from side) 


The main aim of the present paper is to analysis of redistribution of dopant and radiation defects to 
determine conditions, which correspond to decreasing of elements of the considered filter and at the 
same time to increase their density. At the same time we consider a possibility to decrease mismatch- 
induced stress. 


Method of Solutions 

To solve our aim we calculate spatio-temporal distribution of concentration of dopant in the 
considered heterostructure and analyzed it. The above distribution was obtained as solution of the 
following the second Fick's law [1,20-24]. 
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Here C(x,y,z,t) is the spatio-temporal distribution of concentration of dopant; Q is the atomic volume of 
L, 

dopant; V; is the symbol of surficial gradient; |C (x, y, z,t)d z is the surficial concentration of dopant 
0 


on interface between layers of heterostructure (in this situation we assume, that Z-axis is 
perpendicular to interface between layers of heterostructure); 4(x,y,z,t) is the chemical potential due 
to the presence of mismatch-induced stress; D and Ds are the coefficients of volumetric and surficial 
diffusions. The values of dopant diffusion coefficients depends on the properties of materials of 
heterostructure, the speed of heating and cooling of materials during annealing and spatio-temporal 
distribution of concentration of dopant. Dependences of dopant diffusions coefficients on parameters 
could be approximated by the following relations [22-24] 
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Here D, (x,y,z,T) and Dis (x,y,z,T) are the spatial 
(due to accounting all layers of 
heterostruicture) and temperature (due to 
Arrhenius law) dependences of dopant diffusion 
coefficients; T is the temperature of annealing; P 
(xy,z,T) is the limit of solubility of dopant; 
parameter y depends on properties of materials 
and could be integer in the following interval y 
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distribution of concentration of radiation 
vacancies; V" is the equilibrium distribution of 
vacancies. Concentrational dependence of 
dopant diffusion coefficient has been described 
in details in [22]. Spatio-temporal distributions 
of concentration of point radiation defects have 
been determined by solving the following 
system of equations [20,23,24] 
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I (x,y,2,0)=fi (x,y,z), V (xy,Z,0)=fv (x,),Z). (4) 
Here | (x,y,z,t) is the spatio-temporal distribution of concentration of radiation interstitials; I* is the 
equilibrium distribution of interstitials; D)(x,,z,T), Dv(xy,z,T), Dis(xy, 2,T), Dvs(xy,z,T) are the 
coefficients of volumetric and surficial diffusions of interstitials and vacancies, respectively; terms 
V2(x,y,z,t) and I2(x,y,z,t) correspond to generation of divacancies and diinterstitials, respectively (see, 
for example, [24] and appropriate references in this book); kiv(x,y,z,T), ki(x,y,z,T) and kv,v(x,y,z,T) are 
the parameters of recombination of point radiation defects and generation of their complexes. 


Spatio-temporal distributions of divacancies @y (x,y,z,t) and diinterstitials 2 (x,y,z,t) could be 
determined by solving the following system of equations [20,23,24] 
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Here Da(x,y,z,T), Dov(Xy,z,T), Dais (Xy,z,T) and Dovs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of complexes of radiation defects; ki(x,y,z,T) and ky(x,y,z,T) are the parameters of 
decay of complexes of radiation defects. 


Chemical potential z in Eq.(1) could be determine by the following relation [20] 
M=E (Zz) Qo; [Ui(Xy,2,t)+uj(Xy,z,t)]/2, (7) 
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tensor; uj uj are the components u,(x,y,z,t), u,(x,y,z,t) and u,(x,y,z,t) of the displacement vector 
u (x, y, z,t) ; X;, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 
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where ois Poisson coefficient; & = (as-adz_)/az, is the mismatch parameter; a;, az, are lattice distances 
of the substrate and the epitaxial layer; K is the modulus of uniform compression; /is the coefficient of 
thermal expansion; T, is the equilibrium temperature, which coincide (for our case) with room 
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temperature. Components of displacement vector could be obtained by solution of the following 
equations [25] 
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p (Zz) is the density of materials of heterostructure, 6j Is the Kronecker symbol. With account the 
relation for oj last system of equation could be written as 
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Conditions for the system of Eq. (8) could be written in the form 
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We determine spatio-temporal distributions of concentrations of dopant and radiation defects by 
solving the Eqs.(1), (3) and (5) using standard method of averaging of function corrections [26]. 
Previously, we transform the Eqs.(1), (3) and (5) to the following form with account initial 
distributions of the considered concentrations: 
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Farther, we replaced the concentration of dopants and radiation defects in right sides of Equations 
(1a), (3a) and (5a), on their not yet known average values ap. In this situation, we obtain equations 
for the first-order approximations of the required concentrations in the following form 
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Integration on the left and right sides of the Equation (1b), (3b) and (5b) on time gives us possibility to 
obtain relations for above approximation in the final form 
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We determine average values of the first-order approximations of concentrations of dopant and 
radiation defects by the following standard relation [26] 
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Substitution of the relations (1c), (3c) and (5c) into relation (9) gives us possibility to obtain required 
average values in the following form 
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We also determine approximations of the concentrations of dopant and radiation defects 
second and higher orders of concentrations of we replace the required concentrations in the 
dopant and radiation defects framework Eqs. (1c), (3c), (5c) on the following sum anpt+p 
standard iterative procedure of method of ni (xXy,z,t). The replacement leads to the 
averaging of function corrections [26]. following transformation of the appropriate 
Framework this procedure to determine equations 


approximations of the n-th order’ of 


OC,(x,y,zt)_ 8 {regltcralessM ||, Mere), el . 


Ot Ox Piaya.T) Vv" = (vy 


2 
Di (syn2e7) DEN, 0 ae pg Palen 


, Diana tg etele sit}, 0 [1+ Vi yet) Sel 
& 


P’ (x, y,z,T) (v°y 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.10 No.1 3203-3223 (2023) 3212 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Di (x,y, ereale.net)h +é lane + Glu» zot}f + fo (x,y, z)d(t) 


P’ (x, y,z,T) 


L, 


L. 

= L(x, y, Zt a +clsvat law} 

2 |Pev le ML 2C 
0 | D, 

02 | Pevsuls y, Ls I [er rclanwnlaw| 


O1,(x,y, 2,1) O O1,(x, y,z,t) O 1,(x, y,z,t) 
= D T D T }—————— 
At ral 105 95267) By 105 95267) 


6 EsACe ee 
2) Dlssar) fre kv y.zT)a, + h(x, y.z,t)loy +Vi(x y. z,0)] 


Te 
ky HZ T ley +(e y.2.0F +O {2s Vu ( X,Y, Zt ft a, + 1,( x, y,W aqaw| 
0 


0 |D L, 
02 {Pav (x,y, z.t) laa, + h(x, swalaw} (30) 
0 
OV,(x, y, z,t) O OV,(x, y,z,t) Oo OV, (x, y, Zt) 
ate ee a yg a Ed 
Ot ral WaT) a er (x, 9.2.7) = 


e PVE 
2) Dlssar) al k, v(x, y,z,T lay, +(x yz t)llay +Vi y,z,)] 


L, 
ky yl yr2eT lay +Y (9.201) 40 {Be Vt (xy, Zt ft Ay, + V,(x, y.W, slaw 
0 


kT 


2 Dai (X yz) _ £ Ip, (x, y, z, qe) (sy, =] r <I, (x, y, i es us 2 
y 


L, 
a2 {Ps Vu (x, YZ, ‘) [eay + V(x, y,W,t)| aw| 


Ot Ox Ox y 


6 | Dz, 
Uae | kT sH (% Y, Zt flere, +®,,(x, y,W Haw] ek abs y,Z, T)I?(x, y,z,t) 


D ‘f 
vad a V sl (x, yizat){ a0, +O,(s. jaw |ebls mets y, Zt) 
y 0 


2 <| Dg, (xy, er) Ould) + fo, (x y,z)5(t) (5d) 
OZ OZ 


O®,,(x, y,zt) 2 OO, (x, y,zt)| 2 O®,,(x, y, z,t) 
= D .Y,Z51 Dy, (x, y, 2, T ) 
Ot Ox oy (BT) Ox oy (BT) oy 


Pankratov 


International Journal of Advanced Science and Engineering 


www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.10 No.1 3203-3223 (2023) 3213 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


a 
21M vu la ricifltn +OulsrWAaw|hysls et) Vo mcd 


a L, 
oe V sll Gs .208){le, +O,(x,y,W.0) aw stylet V (x,y, z,t) 


+2 1p, (x, y,z,T) on » a] 
v6 


Oz aot om (x, y, z)d(t). 


Integration of the left and the right sides of Equations (1d), (3d) and (5d) give us possibility to obtain 
relations for the required concentrations in the final form 
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Average values of the second-order approximations of required approximations by using the following 
standard relation [26] 
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Next, we determine solutions of Eqs.(10), ie. averaging of function corrections we replace the 
components of displacement vector. To required functions in the right sides of the 
determine the first-order approximations of the equations by their not yet known average values 
considered components framework method of a. The substitution leads to the following result 


pc) taltrtt) pe (2)p(2) ne) 
i Ox 
ol SVs 
ple) ere) 4 (a) p(y Pena) 
t Oy 


, 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.10 No.1 3203-3223 (2023) 3216 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


2 
p(T Melero) (e)p(-) Tree) we 


Integration of the left and the right sides of the above relations on time t leads to the following result 
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Approximations of the second and higher orders of components of displacement vector could be 
determined by using standard replacement of the required components on the following sums 
ajtu;(x,y,Z,t) [26]. The replacement leads to the following result 
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Integration of the left and right sides of the above relations on time t leads to the following result 
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According to the theoretical research, the 
dopant concentration, radiation defect 
concentrations and _ displacement vector 
components are calculated using the second- 
order approximation method of the averaging of 
function corrections. This estimate is usually 
good enough to make a good analysis and get 
very useful results. All results were verified by 
comparison with numerical simulation results. 


DISCUSSION 

In this section, we analyse the kinetics of 
additive and radial defect redistribution during 
annealing and under the influence of mismatch- 
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induced stress. <A_ profile of additive 
concentrations in heterostructures is shown in 
Figures 2 and 3 are for diffusion and ion-type 
doping, respectively. These distributions are 
calculated for the effect of the additive diffusion 
coefficient in the epitaxial process greater than 
the substrate. The figure shows that the 
inhomogeneity of the heterostructure makes it 
possible to increase the compactness of the 
transistor. At the same time, it is seen that the 
contribution distribution of the doped part of 
the epitaxial layer is gradually increasing. The 
compact nature of transistors makes it possible 
to make them fast. 
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L/2 3L/4 L 
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Fig. 2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Increasing of number of curve 
corresponds to increasing of difference between values of dopant diffusion coefficient in layers of 
heterostructure under condition, when value of dopant diffusion coefficient in epitaxial layer is larger, 
than value of dopant diffusion coefficient in substrate 
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Fig. 3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to 
annealing time © = 0.0048(L,2+L,2+L,2)/Do. Curves 2 and 4 corresponds to annealing time © = 
0.0057(L,2+Ly2+L72)/Do. Curves 1 and 2 correspond to homogenous sample. Curves 3 and 4 
corresponds to heterostructure under condition, when value of dopant diffusion coefficient in epitaxial 
layer is larger, than value of dopant diffusion coefficient in substrate 


C(x,®) 


0 L, 


Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time 
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Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time. 
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The second consequence leads to a decrease in distribution cannot be detected. If the annealing 
the local heating of materials during the time is long, the doping concentration 
operation of transistors or a reduction of their distribution is also the same. We recently 
size to the maximum stable local overheating. introduced a method to optimize the annealing 
However, in the framework of such a bipolar time [15,25-32]. Within the framework of this 
transistor fabrication method, the annealing of process, we estimated the true distribution of 
additives and/or electrical defects needs to be additive concentrations with a _ stepwise 
optimized. The reason for this optimization is as function (see Figures 4 and 5). In addition, we 
follows. If the annealing time is short, the determine the effective value of the annealing 
additive does not provide any interaction time by decreasing the mean square error 
between the heterostructure materials. In this below. 
case, a change in the ra concentration 
L,Ly L, 
17 
a ae ae ie (x, Yk @)- w(x,y,z)|dzdydx, (17) 
where y (x,y,z) is the approximation function. dopant achieves appropriate interfaces between 
The dependences of optimal values of annealing materials of heterostructure during annealing of 
time on parameters are presented on Figs. 6 and radiation defects. If dopant did not achieve any 
7 for diffusion and ion types of doping, interfaces during annealing of radiation defects, 
respectively. It should be noted that it is it is practicably to additionally anneal the 
necessary to anneal radiation defects after ion dopant. In that situation, the optimal value of 
implantation. One could find the spreading of additional annealing time of implanted dopant 
concentration of distribution of dopant during is smaller, than annealing time of infused 
annealing. In the ideal case distribution of dopant. 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values of 
dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of dimensionless 
optimal annealing time on value of parameter efor a/L=1/2 and €= y= 0. Curve 3 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and ¢= v= 0. Curve 4 is the 
dependence of dimensionless optimal annealing time on value of parameter yfor a/L=1/2 and e= €=0 
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Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which 
have been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the 
dependence of dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each 
other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of 
dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is the 
dependence of dimensionless optimal annealing time on value of parameter for a/L=1/2 and e= y= 
0. Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for 
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Fig. 8. Normalized dependences of component u; of displacement vector on coordinate z for nonporous 


(curve 1) and porous (curve 2) epitaxial layers 


2000 5 


1000 4 


Wu(x.©) 


0 


0.00 0.25 


T 1 
0.50 0.75 1.00 


Fig. 9. Normalized distributions of charge carrier mobility in the considered heterostructure, the curve 
1 corresponds to the heterostructure, which has been considered in Fig. 1. Curve 2 correspond to a 
homogenous material with averaged parameters of heterostructure from Fig. 1 


Further we analyzed influence of relaxation of 
mechanical stress on distribution of dopant in 
doped areas of heterostructure. Under the 
following condition a< 0O one can find 
compression of distribution of concentration of 
dopant near interface between materials of 
heterostructure. On contrary (at 4>0) one can 
find spreading of distribution of concentration 
of dopant in this area. This changing of 
distribution of concentration of dopant could be 
at least partially compensated by using laser 


annealing [29]. This type of annealing gives us 
possibility to accelerate diffusion of dopant and 
another process in annealed area due to 
inhomogenous distribution of temperature and 
Arrhenius law. Moreover, accounting relaxation 
of mismatch-induced stress in heterostructure 
could leads to changing of optimal values of 
annealing time. Mismatch-induced stress could 
be used to increase density of elements of 
integrated circuits. On the other hand one could 
leads to generation dislocations of the 
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discrepancy. Fig. 8 shows distributions of [8] 

component of displacement vector, which is 

perpendicular to interface between layers of 
heterostructure. 

CONCLUSION 

In this paper, we simulate the implantation and 

redistribution of implanted dopants during the [9] 

fabrication of the field-effect transistors (FET) 

in the framework of a absorptive single pole 

single throw switch, considering the stress 

caused by resting mismatch. We have developed 

recommendations for optimizing annealing to 

reduce the size of transistors and increase their 

speed. We provide recommendations to reduce [10] 

the stress caused by inconsistencies. Analytical 

methods are presented to simulate diffusion and 

ion-type doping, including simultaneous spatial 

and temporal changes. At the same time, the 

method gives us the ability to consider the [11] 

nonlinearity of the decision process. 
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